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It is usually assumed that the boundaries do not affect the bulk microstructure of an 
^ ■ interacting Bose gas. Therefore, the models use the most convenient cyclic boundary 

^ . ' conditions. We show that, in reality, the boundaries affect strongly the bulk microstructure, 

by changing the ground-state energy Eq and the energy of quasiparticles E{k). For the latter, 

we obtain the formula E{k) ^ \/ {^^^Y ~^ ^'2^^ ^^m' differing from the well-known Bogolyubov 
' formula by the factor 1/2^, where f is the number of noncyclic coordinates. The influence 

of boundaries is related to the topology, namely to the difference of the sets of oscillatory 
eigenmodes for closed and open systems. 
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■ 1 Introduction 

^ ■ 
O 

CN ' A weakly interacting Bose gas (WIBG) has been well studied till now (see, e.g., the review 

[1]). The foundations of the microscopic model of WIBG were developed about six decades 
. ago in the classical works by N. Bogolyubov [2] and Bogolyubov and D. Zubarev [3j. Later, 

' the formulas for WIBG were reproduced in numerous works on He II. About twenty years, 

WIBG is studied experimentally with the use of systems containing > 10^ atoms of alkaline 
metals in traps However, we will show that a surprises are possible in this field. In the 
present work, we study the influence of walls of a vessel on the bulk microstructure of WIBG. 

It is commonly considered that the boundaries affect only a microstructure of the near- 
wall layer. Therefore, the convenient cyclic boundary conditions (BCs) are usually used in the 
modeling of bulk properties. But such conditions are realized seldom (and they are impossible 
for the three-dimensional systems, unless a gas occupies the whole Universe). However, though 
the boundaries are far from the majority of atoms in a vessel, their influence on these atoms 
can turn out significant. In fact, we know the "butterfly effect" in complicated nonlinear 
systems, where a small action can transit the system in a qualitatively different state. In work 



[5], the role of boundaries was studied for a Bose fluid by means of the calculation of the wave 
functions (WFs) under zero BCs. It was proved that the boundaries affect strongly the ground- 
state energy of the system, Eq, and the dispersion curve E{k). This result is of importance, and 
it contradicts the commonly accepted expectations. However, in the calculation [5] the WFs 
of the system were expanded in the sets of collective variables p^, which are not independent. 
This is quite admissible, but may cause questions. Therefore, we will solve the problem within 
another method, where WFs are expanded in independent basis functions. 

It was found in several works that the presence of boundaries does not affect the bulk 
microstructure. The main reason for the effect to be missed is the simplified modeling of 
interatomic interaction (e.g., the point interaction). This is discussed in [S', "B]. Possible 
experiments and the general discussion of the influence of boundaries are considered in [5]. 
Below, we consider a system with a nonpoint interaction. 

2 Ground state of a Bose gas 

The boundaries influence a microstructure of the system a) directly by means of the influence 
of the near-surface layer on the bulk properties and b) indirectly due to a difference of the 
collections of harmonics in the expansion of the interatomic potential for closed and open 
systems. It is clear intuitively that the first effect is slight, and the analysis in [Sj confirmed 
this assumption. Therefore, we do not consider the factor in WFs that ensures vanishing of 
WFs on the boundaries. We will consider only the second circumstance, namely the difference 
of the collections of harmonics. 

The key point is simple (see [5] for details). If Xi,X2 G [—Lx/2, Lrc/2], we have Xi — X2 G 
[— L^./2] under cyclic BCs and Xi — X2 G [— L^,La.] under zero BCs. So, under cyclic 
BCs, the expansion of the interatomic potential takes the form 
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where 27r over the sum means that k runs the values 




(3) 



with integers jx,jy,jz- But, under zero BCs, the expansion reads 




(4) 



k 



where tt over the sum indicates that k runs the values 




(5) 



Here, / = 3 is the space dimensionality, and V = L^LyLz is the volume of the system. Since 
the realistic potentials differ from zero significantly only at small |ri — r2|, the values of i>{k) 
are identical for both BCs, though the integration volume over r in ([2]) is different by 2^^ 
times. The difference of expansions (j4j) and ([T]) will cause an effect. 

We note that the potential U{ri — can be expanded not only according to Eq. 
but in different other complete sets of basis functions, e.g., in the functions e«kiri+ik2r2 -^yj^j^ 
ki,k2 ([3]) and where ri and r2 are considered as two independent radius- vectors, or in the 
sines sin kxX sin ky-y sin kzZ, etc. We can use any of these complete collections of functions: 
the expansion in each of the collections reproduces exactly the potential U{ri — r2) under 
the condition that the points riandr2 lie inside the system. In this case, under zero BCs, the 
expansion in the functions e«kiri+jk2r2 jg j^^^ reduced to (jlj) (where ri — r2 is treated as one 
radius- vector). We will use the expansion (jlj), for which the solution can be found by a most 
simple way. 

Let us consider interacting Bose particles, which are placed in a box X Ly X Lz with 
noncyclic boundaries. The Hamiltonian of the system reads 

We seek the ground-state WF in the standard form |T] 

where the prime over the sums means that ji ^ for each rj^ — rj^. In ([7j), the last term 
is the sum with Sn describing the correlations between all particles. Performing the Fourier 
transformation of the functions Sj and extracting the lowest sums from the highest ones, we 
obtain 

kiT^O ■ jij2 
^ ^ a3[ki,k2) giki(r,^-r,2)+ik2(r,,-r,3) ^ 

ki,k2^0 ' iij2i3 

^ ^ a^(ki,0,k3) ^ giki(r,,-r,g+ik3(r,3-r,J ^ 

ki,k3^0 ■ jijzjsji 

k27^ki,k3 /'I, 1, 1, \ ' 
^ ^ Q4lKi,K2,K3j ^ giki(r,,-r,2)+ik2(r,2-r,3)+ik3(r,3-r,J^___^ j^g^ 

ki,k2,k3^0 ' jijzjsji 

where a2(— k) = 02 (k) and 

a3(ki, k2) = a3(-k2, -ki) = a3(-ki + k2, k2) = a3(ki, -k2 + ki). (9) 



In the last term of ([8]) ki 7^ ks. We extract the sum with 04, because the equations yield 
a4(ki,k2 — )■ 0, ka) 7^ a4(ki,0,k3). In expansion ([8]), the lowest sums cannot be extracted 
from the highest ones. Eq. ([HD is the expansion of the logarithm of \I'o in independent basis 
functions, 

1, e'^^^'^h-^'h') ^ g«ki{rj^-rj2)+ik2{rj2-rj3)^ 

g«ki(rjj-rj2)+ik3{rj3-rj4)^ ^iki (r^^ -r^-^ )+ik2 (rj2 -r^g )+jk3 (r^g -rj_j ) ^ ^ ^ ^ j-j^g^ 

(here kj ^ and k2 7^ ki,k3; moreover, ki 7^ ks in e*''^^''^i~'"^2)+«''2(rj2-rj3)+ik3(rj3-rj^)^ 
([8])). Their independence is related to the different numbers of vectors or to the difference 
in ji or in k^ (for the same r^). The independence holds under cyclic and zero BCs, since 
the exponential functions are the basis functions of a Fourier expansion. It is of importance 
that the wave vectors k; enter fllOl) at the difference of coordinates. Therefore, k; run values 
([5]) and ([3]) under zero and cyclic BC, respectively. Since the functions contain the k- and 
r- variables, such an approach can be called the "fcr-method". Similar approaches were used 
earlier: structure ([7]) was considered by Feenberg [7], Vakarchuk and Yukhnovskii and 
Krotscheck [9]. However, expansion ([8]) (as well as ( 12T|) and (122|) below) was not applied: 
either another representation was used [7J; or the terms with = ji were introduced ([8]); or 
the transition to the collective variables pk was made [8J; or the study was carried out only 
in the r-representation [H]. We proceed from structure ([HD, (1211) . and (I22p . where the basis 
functions are independent for various BCs. Since the coefficients in each sum in ([8]) are 
identical for different and identical k^, we will collect the coefficients of sums 

(which are the collective variables), rather than of separate functions ffTOj) . 

We now pass to the solution. For \E'o (H) and Hamiltonian (jS]), the Schrodinger equation 

= iV^o^o (12) 

is reduced to 

- ^ + (^^•^)'] + ^ £ - r.l) = A^^o. (13) 



Let us substitute (jlj) and ([8]) in the last equation. In view of the independence of functions 
( JTOj) . Eq. (fT3|) holds, if the coefficients in front of sums ( ITT]) and of the constant are equal to 
zero. After calculations, we obtain the equations for the ground-state energy Eq (per atom) 



and the functions af 



2/+1 2m ' N ^ 2m 

k^O ki,k27^0 



a2(k)e-al(k)e = i-^a2(q){(g' + kq)a2(k + q) + 



q^O 

+ (2g2_kq)a3(k,q) + (g2 + kq)a3(q,k + q) + g2a3(k + q,q) } + ..., (15) 

a3(ki,k2)[eo(ki) + eo(k2) + eo(k2 - ki)] = 2kik2a2(ki)a2(k2) + 
+ 2k2(k2 - ki)a2(k2)a2(k2 - ki) + 2ki(ki - k2)a2(ki)a2(ki - k2) + 
6 

+ ^ $^(fc3 + k3k2)a2(k3)[a3(ki,k2 + k3) + a3(ki + k3,k2 + k3)] + ..., (16) 

k3^0 

6o(k) = k2(l-2a2(k)), (17) 

where the dots mean the higher corrections (the sums with a1, a^, a% in ffl^ and the sums 
with Og, 02*^4, 0'20'4:, etc. in (1151) and (fTB]) ). We omit the equations for the functions aj>4. 

For a weak interaction, the corrections with aj>3 can be neglected. With regard for (fT9|) . 
relations ([I1D-([I7]) yield 

^ nz/(0) 1 p. 1 p. nvik) , , 

'^^ + a,{k)k^-alik)k^ = 0. (19) 

The same formulas were obtained in [5j within a different method. 

For cyclic BCs, is is necessary to use expansion ([1]). This leads to the replacements 

(tt) (27r) 

v{k)/2f —7- v{k), XI X] equations, and formulas (lT8l)-( IT9|) become the well-known 

Bogolyubov formulas. Formula ( IT^ for cyclic BCs without corrections with aj>z was deduced 
earlier in [7]. 

3 Excited state with one quasiparticle 

The WF of a state with one phonon is sought in the form |10[ |3] 

^k(ri,...,r^) =V'k^o, (20) 
= P-k + ^V'k, (21) 

N 

where pk = X] ) and ^^/^k corresponds to higher corrections. Their exact form is 
determined by the Schrodinger equation. We write them similarly to (jH]): 

{^) 

JV-k = &2(qi,k)5^e'('^+'ii)'-^i-*'ii'-^2 + 

qi7^o,-k jij2 

(7r)qi+q2+k7^0 / 

+ 63(qi,q2,k) e^^'^^'"'^'''^'''"*'''''' + ••• (22) 



Here, 



&2(q,k) = 62(-k-q, k). 



(23) 



h{<ii, q2, k) = 63(q2, qi, k) = &3(-qi - q2 - k, q2, k) = b^iq^, -Qi - q2 - k, k). (24) 

This relation differs from the solution in the p^-representation [TT| |5] by the absence of terms 
with jk = jl in sums. Therefore, the higher sums cannot be reduced to the lower ones (it 
is possible in the pk-representation, since the terms with = ji have a structure of lower 
sums) . 

The Schrodinger equation yields the equation for ip]^: 

-^Yl {^^-^k + 2(V,-^k) X V,S} = E{k)^^, (25) 

j 

where -E'(k) is the phonon energy. Let us substitute ip^ (12T]) . ( !22|) in ( !25|) . Formulas ( 12T|) and 
(122]) give the expansion of ip^ in the independent (for a fixed k) functions 

(their independence can be easily seen if we take into account that, for example, 

gi(k+qi+q2)rj;^-iqirj2-iq2r^3 ^ gikr^-^ +tqi (r^-^ -rj^ )+jq2(rj;^ -r^g ) ^ giq(rj^-rj2) g^j.g ^]^g basis fuUCtioUS 

of the expansion of a function of the form /(r^^ — tj^) in the Fourier series). Equation (l25l) 
holds, if the coefficients in front of the functions (l26l) are zero. This yields the equations for 
E{k) and bji 

e(k) = eo(-k)-^ &2(q,k)[g2a2(q) + (g2 + kq)a2(_q)] _ 

q^O,-k 

- ^ 5Z ?'^2(q,k)a3(q,-k)-^5^g'a2(q)fc3(q,-k,k) + ..., (27) 

QT^O-k q^O 

62(q, k) [eo(q) + eo(-k - q) - e(k)] = -kqa2(q) + k(q + k)a2(-q - k) + 
+ A;2a3(q + k,k)/2 + A;2a3(-q,k)/2 + |^ X (28) 
X J2 a2(qi){(gi+qiq)fe2(qi + q,k) + (gi2-qi(q + k))62(qi-q-k,k)} + ..., 

qi7^0,-q,-k-q 

&3(qi,q2,k)[eo(qi) + eo(q2) + eo(-k-qi - q2) - e(k)] ^ 
~ -k(qi + q2)a3(qi + q2, q2) + k'^a^ik + qi + q2, k + q^, k)/3 + 
+ 2(qi + q2)62(qi + q2,k)[2qia2(qi) + (qi + q2)a3(qi + q2,q2)]- (29) 

Here, e(k) = 2m-E(k)//i^, and the dots mean the higher corrections (the sums with ai±ibi and 
aihi in (127|) and the sums with ai±2bi, ai±ibi, and aibi in (128|) ). We omit the equations for bj>4 
and write the equation for 63 in the zero approximation (we do not symmetrize its right-hand 



side according to ( IMl) to avoid too bulky formulas). We also neglect in coefficients of the 
form 1 + in Eqs. (IH-([I7D and i^-l^. 

In the presence of walls, the stationary state is a standing wave, rather than a traveling 
one. The solution for a standing wave is a superposition of traveling waves (the normalizing 
factor is omitted), 

'i/'k = V'k + 7 permutations, (30) 

where the permutation means ip^ with the different sign of one or several components of the 
vector k. The energy of such a wave coincides with that of a traveling wave. If the interaction 
is switched-off, \l/k must be reduced to the solution for a free particle in box (in this case, it is 
necessary to consider the vanishing of the WF on boundaries [5J). For this solution, k takes 
values ([5]). Therefore, k of a phonon must take also the same values. 

For a weak interaction, we can neglect the corrections with bj. Then we have 

Taking (I19p into account, we obtain the Bogolyubov formula 



but with the additional factor due to boundaries. This factor is absent under cyclic 
boundaries. In the general case, / is the number of noncyclic coordinates. Formula ( l32l) was 
also obtained in [5] within the pk-niethod. 

Note that it is easy to see without any calculations that the transition from cyclic to zero 
BCs requires one to add the factor 1/2-^ at z/(fc) in the formulas for Eq and E{k): indeed, the 
potential enters the Schrodinger equation only as sum (jlj), in which the factor 1/2'^ arises. 

In the proposed fcr-method, each of the equations in the chains of equations for \l/o and 
ipi^ contains ~ A^^ terms (but only several terms in pk-method [TT| E]). It is impossible 
to solve such equations for a strong interaction. However, the fcr-method is effective for a 
weak interaction. The advantage of the A;r-method consists in the independence of the basis 
functions not only under cyclic BCs, but under zero BCs as well. 

Since \l/o is quantitatively different for zero and cyclic BCs, the amounts of a condensate 
must differ significantly for these BCs, especially for Bose systems with high density. 

In [6j, the problem with boundaries (ID, / = 1) was solved within the Gross-Pitaevskii 
approach, and two solutions for the dispersion law were found: traditional Bogolyubov solution 
and the above-obtained one (132|) . Two solutions mean the existence of two ways to the 
diagonalization of a Hamiltonian, with two different ground states. But, in the approach 
with A^-particle WFs, two ground states are forbidden by the theorem of nodes (the case 
with degeneration is improbable). One solution for the ground state was found above and in 
[5]; it corresponds to the dispersion law ( !32|) . Hence, the "ground state" for the traditional 



dispersion law |6] has nodes and is not the lowest state of the system. In this case, the 
traditional dispersion law must not be realized, if / > 0. 

4 Conclusion 

We have shown that the boundaries affect strongly the bulk microstructure of a Bose gas 
and a Bose fluid. This influence is not quite trivial and is related to the difference of the 
collections of oscillatory eigenmodes for closed and open systems. As far as we understand, 
the effect is due to the fact that a standing wave keeps the memory of walls, because it is 
modulated by them. But a particle loses the memory of a wall after several collisions with other 
particles. Therefore, this effect must be typical of all weakly excited systems, whose physics 
is determined by collective excitations (waves), rather than separate atoms. In particular, it 
must be inherent in quantum crystals, superconductors and some Fermi systems, for which 
the Bogolyubov dispersion law was obtained |12] . 
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